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We carried out first-principles calculations within Density Functional Theory to investigate the 
structural, electronic and magnetic properties of boron-nitride (BN) honeycomb structure function- 
alized by adatom adsorption, as well as by the substitution of foreign atoms for B and N atoms. 
For periodic high density coverage, most of 3d transition metal atoms and some of group 3A, 4A, 
and 6A elements are adsorbed with significant binding energy and modify the electronic structure 
of bare BN monolayer. While bare BN monolayer is nonmagnetic, wide band gap semiconductor, 
at high coverage of specific adatoms it can achieve magnetic metallic, even half-metallic ground 
states. At low coverage, the bands associated with adsorbed atoms are flat and the band structure 
of parent BN is not affected significantly. Therefore, adatoms and substitution of foreign atoms at 
low coverage are taken to be the representative of impurity atoms yielding localized states in the 
band gap and resonance states in the band continua. Notably, the substitution of C for B and N 
yield donor and acceptor like magnetic states in the band gap. Localized impurity states occurring 
in the gap give rise to interesting properties for electronic and optical application of the single layer 
BN honeycomb structure. 

PACS numbers: 73.22.-f, 73.90.+f, 75.50.Pp, 75.75.-c 



I. INTRODUCTION 

Research on BN based materials have grown gradually 
in recent years P"^ This is not only due to their fascinat- 
ing properties, such as hardness, high melting point, and 
large band gap, but also due to the geometric similarity of 
planar, two dimensional (2D) BN to graphene. Scientists 
already achieved the synthesis of single layer BN hon- 
eycomb structure on substrates^' and a few layer thick 
structures from 3D hexagonal (h-)BN either on a sub- 
strate or freestanding PES' Recently, Jin et a/.Hl reported 
the fabrication of freestanding BN honeycomb structure 
(we specify it simply as 2D BN throughout the text). The 
realization of the synthesis of 2D BN is rapidly attracting 
interest on BN, since it has 2D hexagonal lattice, which 
is commensurate to the lattice structure of covalently 
bonded graphene. More recently, the synthesis of sin- 
gle layer composite structures consisting of adjacent 2D 
BN and graphene domains is realized.^ However, unlike 
semi- metallic graphene, 2D BN is a nonmagnetic, wide 
band gap semiconductor with an indirect energy gap of 
4.64 eV^ calculated within generalized gradient approx- 
imation (GGA). The indirect gap is further corrected to 
6.82 eV with GWo self-energy method by §ahin et al™A 
theoretical comparative study of 3D and 2D BN, and its 
nanoribbons comprising their mechanical, electronic and 
magnetic properties was reported by Topsakal et a/P^ 
2D BN, and their nanoribbons can be easily functional- 
ized by mairy different ways for different purpo ses su ch 
as doping^tS, exchange of atoms and vacancies! 18 ' 19 * 

In this paper, using state-of-the-art first principles 
plane wave calculations we investigate the effects of 
adatoms adsorbed on 2D BN, as well as the substitution 
of foreign atoms for B and N atoms in the honeycomb 
structure. We consider both high coverage (where the 



coupling between adjacent foreign atoms is substantial) 
and low coverage (where the coupling is negligible). We 
conclude that at high coverage (or decoration) of spe- 
cific adatoms one can turn the nonmagnetic, wide band 
gap material into magnetic, metallic or even half- metallic 
states. At low coverage, adatoms give rise to localized 
states in the band gap. 



II. METHOD 

We perform first-principles, spin-polarized plane 
wave calculation d 20 ' 21 ' within density functional the- 
ory (DFT) 22 using PAW (projector augmented- wave) 
potentials. 23 The exchange-correlation potential is ap- 
proximated by generalized gradient approximation 
(GGA)! 2 -^ For the partial occupancies, we use the 
Methfessel-Paxton smearing method!^ The width of the 
smearing for all structures is chosen as 0.01 eV for ge- 
ometry relaxations and band structure calculations. For 
accurate density of states (DOS) calculations the width of 
smearing is taken as 0.1 eV. We consider a single adatom 
adsorbed to each (2x2) and (4x4) supercells of 2D BN 
structure and treat the system using periodic boundary 
conditions. For high coverage corresponding to O = 1/8, 
(2x2) supercell is used , while low coverage, O = 1/32, is 
treated by using (4x4) supercell. A large spacing (at least 
~ 14 A) between adjacent 2D BN layers is taken to pre- 
vent interlayer interactions. The number of plane waves 
used in expanding Bloch functions and k-points used in 
sampling the Brillouin zone (BZ) are determined by a 
series of convergence tests. In the self-consistent poten- 
tial and the total energy calculations, the BZ is sampled 
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Figure 1. (Color online) A (4x4) supercell of 2D BN honey- 
comb structure, which consists of four (2x2) supercells. Var- 
ious possible adsorption sites of adatoms are indicated. The 
hollow site, H: the adatom (grey/orange ball) is placed on 
top of the center of a hexagon. TN site: the adatom is placed 
on top of nitrogen (dark grey/red). TB site: the adatom is 
placed on top of boron (light grey/green). The bridge site BR: 
the adatom is located on top of the boron-nitrogen bond. 



by (15x15x1) mesh points in k-space within Monkhorst- 
Pack scheme^ for the (2x2) supercells. For calcula- 
tions involving (4x4) supercells, the number of k-points 
is taken as (9x9x1). For accurate density of states (DOS) 
calculations, k-points samplings are further increased to 
(25x25x1) and (15x15x1) for the (2x2) and (4x4) super- 
cells, respectively. A plane- wave basis set with the ki- 
netic energy cutoff h 2 \k + G| 2 /2ra = 520 eV is used. All 
the atomic positions and lattice constants are optimized 
by using the conjugate gradient method, where the total 
energy and the atomic forces are minimized. The conver- 
gence is achieved when the difference of the total energies 
of last two consecutive steps is less than 10 -5 eV and the 
maximum force allowed on each atom is 0.03 eV/A. The 
pressure on the system is kept smaller than ~ 1 kBar 
per unit cell in all of the calculations. As a rule, the 
structure becomes more energetic as its total energy is 
lowered. Charge transfer values are calculated according 
to the Bader analysis.^ 



III. ADSORPTION 

The lowest energy sites of various adsorbed atoms are 
determined by placing foreign atoms initially to four pos- 
sible adsorption sites at a height of ~ 2 A from BN plane 
as described in Fig. [T] Upon fully self-consistent ge- 



ometry optimizations with both spin-polarized and spin- 
unpolarized configurations, where all atoms in the super- 
cell are relaxed in all directions, we determine the equi- 
librium site as the lowest energy configuration among 
four different sites. The binding energy of an adsorbed 
adatom is defined as E^ = Ebn + Ea — Ebn+a, where 
Ebn is the total energy of bare 2D BN, Ea is the total 
energy of free adatom calculated in the vacuum. Ebn+a 
is the total energy of 2D BN structure with adsorbed 
adatom. We investigated the adsorption of following sin- 
gle adatoms, namely Sc, Ti, V, Cr, Mn, Fe, Mo, W, Pt, 
H, C, Si, B, N, O, Ca, Cu, Pd, Ni and Zn. Among these 
atoms, Cr, Mn, Mo, W, H, N, Ca and Zn cannot bind to 
2D BN. 



A. Adsorption of Adatoms to (2x2) BN supercell 

(e = i/8) 

The electronic and magnetic properties of 2D BN are 
modified through adatom adsorption at high coverage. 
The adatom-adatom distance is ~ 5 A at O = 1/8, where 
the magnetic coupling may be crucial in determining the 
ground state. In order to account for the antiferromag- 
netic (AFM) coupling between adatoms and to allow 
their reconstruction we treated 8 = 1/8 coverage in (4x4) 
supercell, which includes four (2x2) supercells each hav- 
ing a single adatom. This way adatom-adatom distance 
of ~ 5 A is maintained. We further carried out geometry 
relaxation with three different initial magnetic ordering 
of adatoms. First case is where adatoms are coupled anti- 
ferromagnetically. The second case is similar to first one, 
but adatoms are initially coupled ferromagnetically(FM). 
The last case corresponds to a spin-unpolarized, nonmag- 
netic (NM), geometry relaxation. Our calculations indi- 
cate that Cu, Ni, Pd and Pt adatoms have nonmagnetic 
ground state for O = 1/8, whereas Sc, Ti and V have 
AFM ground state. Oxygen is the only adatom which is 
found to be in the FM ground state. Boron and carbon 
adatoms undergo a reconstruction to lower the total en- 
ergy in (4x4) supercell at O = 1/8. Among all adatoms, 
Si and Fe present the most interesting situation, where 
2D BN monolayer is changed to half- metallic state. How- 
ever, these adatoms have small binding energies and thus 
are excluded from our further analysis. Also Cu having a 
binding energy smaller than 0.25 eV is also excluded from 
our study. Our results are given in Table [T| In Figure [2j 
the calculated energy band structures and correspond- 
ing total density of states (TDOS) and partial density of 
states (PDOS) of Ni, Pd or Pt adatom+2D BN system 
are presented. The band structure of bare 2D BN folded 
to (2x2) BZ is also presented to reveal the effect of the 
adatom adsorption on the electronic structure. 

Nickel atom with As 2 + 3d 8 electronic configuration is 
adsorbed at TN site. Even if Ni atom is placed at H 
or BR sites, it eventually moves to TN site. Flat bands 
slightly below the Fermi energy are all derived from the 
3d orbitals of Ni atom. The band above the Fermi level 
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Table I. Calculated structural, electronic and magnetic properties of 2D BN monolayer at uniform = 1/8 adatom coverage. 
Equilibrium positions of the adatoms, such as TN, TB, H and BR are described in Fig. [l] distances of adatom to the nearest 
N, (In (in A); distances of adatom to the nearest B, ds (in A); the height of adatom from the BN plane, h (in A); the average 
B-N bond length, dsN (in A); the binding energy of adatom, Eh (in eV)\ the net magnetic moment per supercell, p (in /jlb). 
Electronic structure is specified as metallic (M) or semiconductor (SC). The type of the band gap can be either direct (dr) or 
indirect gap (id). The energy gap of the system after adsorption, E g (in eV); the transfer of charge to the adatom from BN, Ap 
(in electrons and Ap < 0, if adatom is negatively charged); the dipole moment of the system along the z— direction, p (in eA); 
work function (or photoelectric threshold for semiconductors) is <I> (in eV). For the adatoms, which have either nonmagnetic 
(NM) or ferromagnetic (FM) ground states, calculations are performed using (2x2) supercell; whereas for adatoms, which have 
either antiferromagnetic (AFM) ground state or reconstruction, (4x4) supercell are used with = 1/8. For adatoms undergoing 
a reconstruction, Ap is given for the average charge transfer per adatom. 
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is mainly formed from the combination of 4s and 3d z 2 
orbitals of Ni. Accordingly, the band structure of parent 
BN layer is not distorted considerably. The charge trans- 
fer between Ni and 2D BN is relatively small and does 
not cause any significant dipole moment. 



Adsorption of Pd and Pt give rise to electronic and 
magnetic properties similar to those of adsorbed Ni. 
They are adsorbed also at TN site and have high binding 
energies. Since the radius of Pd and Pt are relatively 
larger than Ni, the distance to the nearest N atom of 2D 
BN, dj\[ is slightly larger. The charge transfer between 
2D BN and the adatom is small. Ni, Pd and Pt have ion- 
ization energies of 7.63, 8.33 and 8.96 eV, respectively.^ 
This ordering of the ionization energies complies with the 
ordering of work function, <£, of Ni, Pd and Pt covered 
2D BN. For both Pd and Pt, the bands slightly below the 
Fermi level are all derived from localized d orbitals, but 
the band above Fermi level is mainly formed from s or- 
bital of adatom with some d z i contribution similar with 
the case in Ni adsorption. Even if Pd and Pt have the 
same indirect band gap energy, Pd shows more dispersive 
band slightly above the Fermi level. 



TDOS and PDOS corresponding to © = 1/8 cover- 
age of B, C (which undergo a reconstruction in (4x4) 
supercell); Sc, Ti, V (which have AFM ground state), 
and O (which has FM ground state) are presented in 
Fig. [3j Four carbon adatoms treated in a (4x4) supercell 



undergo a reconstruction^ while three C atoms are ad- 
sorbed near TN sites, the remaining one is moved to the 
BR site. Upon reconstruction, the charge on C adatom at 
the bridge cite is minute positive, whereas those at TN 
site are negatively charged by ~ 0.6. Carbon adatoms 
at TN site take their charge mostly from the nearest N 
atom. The bonding of C at TN cite has an ionic character 
resulting a dipole moment on the system. 

Sc, Ti and V have AFM ground state at 8 = 1/8P 
Unlike Sc and Ti, adsorption of V does not make the 
system metallic. The energy bands below the Fermi level 
is mainly from d xyi d yzi d xz and above the Fermi level is 
mainly from d x i and d z i. The energy band gap of 0.19 
eV originates from the splitting of 3d orbital states. 

Among all adatoms in Table [TJ only O has FM ground 
state and is absorbed at TB site. The excess charge on 
the adsorbed O is transferred from three nearest N atoms, 
which leads to an ionic character in bonding. Since un- 
derlying B atom is positively charged and O is negatively 
charged, adsorbed O pulls the nearest B atom and causes 
distortion on planar 2D BN monolayer. 

B. Adsorption of Adatoms to (4x4) BN supercell 

(6 = 1/32) 

We next investigate the adsorption of a single adatom 
to (4x4) supercell, which corresponds to low density cov- 
erage, = 1/32. As indicated in Fig. [I] we initially 
placed the adatoms at four different positions and relaxed 
their geometric structure with both spin-polarized and 
spin-unpolarized calculations. Interestingly, we observe 
that among 9 different adatoms, the adsorption sites of 
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Figure 2. (Color online) (a) The energy band structure of bare 2D BN folded to the (2x2) supercell and corresponding total 
density of states (TDOS). Zero of band energy is set at the Fermi energy, Ep- (b)-(d) The energy band structure of single Ni, 
Pd, and Pt adsorbed to each (2x2) supercell of 2D BN (G = 1/8). Corresponding TDOS (continuous black/blue line) and partial 
density of states (PDOS) projected to the adatom (shaded dark/red) are also indicated. The energy gap of semiconductors are 
shaded (light /yellow). The ground states are nonmagnetic. 



B, C and O differ by going from 6 = 1/8 to 6 = 1/32. 
Since the adatom-adatom distance between adjacent su- 
percells is ~ 10 A, it can normally be contemplated that 
the coupling between adatoms are negligible. To verify 
this assumption, we carried out binding energy calcula- 
tions for C, O, Ti, Sc and Pt adatoms in (8x8) supercell 
corresponding to a coverage of 6 = 1/128. In Table [Ti] 
we compared the equilibrium binding sites and binding 
energies, and magnetic states of these atoms for three 
different coverage 6 =1/8, 1/32 and 1/128. 

The binding energies as well as the magnetic ground 
states of C, O, Sc and Ti vary substantially by going 
from 6 = 1/8 to 6 = 1/32. For C and O, the adsorption 
sites are also changed. However the situation is not the 
same when the coverage is further lowered from = 1/32 
to © = 1/128, hence when the adatom-adatom distance 
is increased from ~ 10 A to ~ 20 A. For adatoms in- 
cluded in Tab. [TTJ the binding energies and magnetic mo- 
ments are not changing significantly, and the adsorption 
sites are remaining the same by going from O = 1/32 to 
O = 1/128. This finding is corroborating our arguments 
made at the beginning of the paper, that the properties 
calculated for 6 = 1/32 coverage (or the adsorption of 
a single adatom adsorbed to each (4x4) supercell) can 
mimic the adsorption of single isolated adatom on a very 
large area of 2D BN (or very large adatom-adatom dis- 
tance). Consequently, the bands of adatom+2D BN cal- 



Table II. Binding energies (i?b), magnetic moments O), ad- 
sorption sites and average adatom-adatom distances^ (gU- a) 
of C, O, Sc, Ti and Pt adatoms adsorbed on 2D BN at differ- 
ent coverages. 
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1/8 E » W 
M Ob) 
dA-A =~ 5 A Site 


1.38 0.98 0.90 0.92 1.48 
2.00 0.00 0.00 - 
BR+TN TB H H TN 


1/32 Eb W 

m Ob) 

d A -A =~ 10 A Site 


1.17 2.01 0.43 0.74 1.52 
2.00 - 3.00 4.00 - 
BR BR H H TN 


1/128 Eb < eV > 

m Ob) 

d A -A =~ 20 A Site 


1.15 1.98 0.37 0.70 1.47 
2.00 - 3.00 4.00 - 
BR BR H H TN 



culated at B = 1/32 become rather flat and can be taken 
as the localized impurity state (or resonances if 2D BN 
states are significantly contributed). Under these circum- 
stances, the band gap and the edges of valence and con- 
duction bands can be unaltered. In Table [TTI| we include 
the calculated geometric, electronic and magnetic prop- 
erties of a single adatom adsorbed to 2D BN at = 1/32. 
The energies of localized and resonance states relative to 
2D BN's valence band edge are also tabulated. 



Density of States Density of States 



Density of States 



Figure 3. (Color online) TDOS (dark/blue line) and PDOS 
(shaded grey /red) for B, C, O, Sc, Ti and V adatoms at high 
coverage, = 1/8. Up (down) arrow on the right (left) site of 
TDOS and PDOS indicates spin direction. Zero of the DOS 
energy is set to the Fermi energy, Ef of the adatom+2D BN 
system which is indicated with dashed-dotted grey /green line. 
All adatoms treated in (4x4) supercell to allow antiferromag- 
netic interaction or reconstruction, whereby single adatom is 
adsorbed to each (2x2) supercells of (4x4) supercell amount- 
ing O = 1/8. Sc, Ti and V have antiferromagnetic ground 
state; O has ferromagnetic state; B and C are nonmagnetic. 
The state densities are given in arbitrary units, but have the 
same scale for all adatoms. 



Boron adsorbed on 2D BN on (4x4) supercell exhibits 
electronic, magnetic and structural properties, which are 
significantly different from those of 6 = 1/8. For exam- 
ple, B atom of the planar BN rises ~ 0.2 A; this slightly 
changes the sp 2 -hybridization locally. Two spin-up bands 
and one spin-down band for spin-down are filled below 
the Fermi energy. In Figure [4j the band structure of 
B+2D BN and isosurfaces of charge density correspond- 
ing to the localized states Ei, E3 and E5 are shown. 
The magnetic properties are different from the case at 
9 = 1/8, since reconstruction of B adatoms does not 
take place due to the absence of adatom-adatom interac- 
tion. 

Adsorption of C on (4x4) supercell of 2D BN causes 
significant splitting of degenerate p orbital levels. Spin- 
up and spin-down bands indicated in Fig. [4] originate 
from p x (Ei state) and p y (E2 state) orbitals of adsorbed 
C atom. In this low coverage, C adatom is adsorbed to 
BR position, but closer to N atom. Similar to B, the 
magnetic properties of the C+BN system changed from 
NM to FM upon lowering the coverage density. 



At 9 = 1/8, O adatom creates a surface distortion and 
pulls underlying B atom upwards. However at 9 = 1/32, 
O adsorbed at BR site does not generate a distorted BN 



region. The bands near the edge of valence band have 
significant dispersion due to strong coupling between O 
and 2D BN. The AFM ground states of both Sc, Ti and 
V adsorbed 2D BN at 9 = 1/8 change to FM ground 
states at 9 = 1/32. Among all adatoms studied in the 
present paper, electronic and magnetic properties of Ni, 
Pd and Pt are not affected upon lowering the coverage 
density. 

Next we address the question whether adsorbed 
adatom migrates on 2D BN honeycomb structure. Here 
we consider only single C and Ti atoms as prototypes 
and calculate the energy barrier for their motion along 
the symmetry directions of 2D BN. In Figure |5j the vari- 
ation of energy of single C and Ti adatoms along sym- 
metry directions are presented. The potential barrier on 
the migration path for C adatom is ~0.30 eV. In the case 
of Ti the energy barrier on the possible migration path 
is only 0.15 eV. The latter barrier is not high enough 
to prevent Ti atoms from cluster formation at elevated 
temperatures. 



IV. SUBSTITUTION (9 = 1/32) 

We, finally deal with the substitutional doping of the 
specific atoms. Here we consider that Al, Be and C sub- 
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Table III. Calculated values for single adatom adsorbed to each (4x4) supercell, corresponding to the coverage G = 1/32. 
Adsorption site, binding energy, Eb\ magnetic moment per (4x4) supercell, /x; the distance from 2D BN monolayer, h\ the 
distance from the nearest N atom, cIn] the distance from the nearest B atom, ds] charge transfer, Ap; and energies of relevant 
localized or resonance states measured from the top of the valence band, E n . (Ft) indicates that the corresponding spin-up 
state is full. (E^) represents an unoccupied spin-down state. If no spin direction is indicated, that state is nonmagnetic. R 
is resonance state having significant contribution from 2D BN states in the band continua. Adsorption site of the adatoms 
corresponding to their lowest total energy are indicated by TN, TB, BR, or H as described in Fig[l] 
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stitution for B atom; C, O, and P substitution for N 
atom of 2D BN honeycomb structure. We mimic the 
substitutional doping by a model, where a specific atom 
substitutes for a single B or N in every (4x4) supercell 
of 2D BN. In this model the distance between impurity 
atoms is ~ 10 A resulting in negligible coupling between 
them. The substitution energy E s is calculate d 32 1 33 1 as: 



E* = 



N - 1 



TV 



-Ebn + Ea — E B \j n — Ebns 



Here TV is the number of atoms in the supercell, which 
is N — 32 for (4x4) supercell. Ebn is the calculated to- 
tal energy of 2D BN of the (4x4) supercell corresponding 
to 16 B-N atom pairs. Ea is the experimental cohesive 
ener g^|28] Q £ impurity atom in its equilibrium crystal. 



Ebns is the total energy of the supercell after substi- 
tution process. E B \j n is the cohesive energy of either 
B or N that depends on which atom is exchanged. The 
experimental cohesive energy of B is defined with respect 
to the B crystal! 2 ^ However, in the case of N, the cohe- 
sive energy of N is calculated relative to N2 in the gas 
phase. According to this expression, a positive energy 
value means the exchange of foreign adatom either with 
B or N is endothermic reaction, whereas the negative 
value indicates an exothermic process. Structural relax- 
ations and the lowest energy states of the final structures 
are further tested with three different initial magnetic 
moment distributions on the atoms and all these cases 



are converged to the same values given in Tab. [IV] Here 
the substitution of C atom for B (N) is of particular in- 
terest, since it dopes the 2D BN honeycomb structure as 
donor (acceptor). 



We now examine the electronic and magnetic structure 
of 2D BN substituted by Al, C, and P described schemat- 
ically in Fig. [6] Al substituting for B and being in the 
same group with B, but having a relatively large atomic 
radius, distorts the planar structure of BN layer and is 
located at a position 0.52 A higher than the substituted 
B atom. The valence band of the BN layer does not in- 
fluence much upon substitution, but states derived from 
p-orbitals of Al adatom appear near the conduction band 
edge. In Figure [6j p^-orbitals of Al are dominant in Ei 



state, however p y is the most contributed in E 2 state. In 
the case of C substituting B, the valence and conduction 
bands of BN layer are not influenced, except minute split- 
ting between the spin-up and spin-down states resulting 
from the magnetic C atom. The excess electron of the 
substituted C atom relative to B leads spin-polarization 
and fills the Ei state in Fig. [6] like a n-type semicon- 
ductor. When substitute for N, C atom creates electron 
deficiency; states derived from p-orbitals of C atom are 
located near the edge of the valence band like a p-type 
semiconductor. p z orbital of C splits and spin-up (E3 
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Figure 4. (Color online) Schematic diagram of the relevant energy levels (or bands) of adatom (B, C, O, Ti and Pt) adsorbed 
to (4x4) supercell of the single layer BN. The light/grey shaded region in the background is the valence and conduction band 
region of 2D BN. Zero of band energy is set at the Fermi level of the parent, bare 2D BN, Ef- Spin-up and spin-down bands 
are shown by dark/blue and light/green lines, respectively. Solid bands indicates that the contribution of adatom to the band 
is more than 50% except O adatom. Energies of some of the relevant adatom (impurity) states relative to valence band edge 
of parent 2D BN are indicated. Below each energy diagram of adatoms, charge density isosurfaces of specific states are shown. 
The isosurface value is taken as 7xl0 -5 electrons/ A 3 . (Note that all localized states of O have significant contributions from 
BN states and are shown by solid dark/blue lines.) 



Table IV. Substitution energy, E s \ magnetic moment per (4x4) supercell, /x; the height of the substituting atom from 2D BN 
monolayer, h; the distance from the nearest B or N atom, d B \j n] charge transferred to substituting atonr^, Ap; and energies 
of relevant localized states in the band gap, E n . (Ft) indicates that the corresponding spin-up state is full. (E^) indicates 
unoccupied spin-down state. If no spin direction is indicated, that state is nonmagnetic. 
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Figure 5. (Color Online) Variation of energy of a single C 
and Ti adatom as they are moving along special directions 
of a hexagon of 2D BN honeycomb structure. The possible 
migration path of adatom facing minimum energy barrier is 
shown by stars on the hexagon. Calculations are carried out 
in (4x4) supercell. 



state in Fig. |6| becomes occupied. Because of this un- 
paired electron C atom attains net total magnetic 
moment in both n-type and p-type substitutional doping. 



Substituted P having relatively larger atomic radius, 
gives rise to the local deformation and raises 1.42 
A above the plane of 2D BN honeycomb structure. 
This gives rise to a local dehybridization of planar sp 2 - 
bonding. Since the upper valence bands of parent 2D BN 
are derived mainly from nitrogen atoms, these bands are 
affected upon their exchange with foreign atoms. Hence 
several resonances appear, such as Ei and E 2 in the up- 
per part of valence band when P substitutes for N atom 
in Fig. [6] In addition a localized state E3 occurs near 
the conduction band edge. 



V. DISCUSSION AND CONCLUSIONS 

In this paper, we showed that 2D BN can be func- 
tionalized to attain properties, which can be useful in 
future applications in nanoelectronics and nanomagnet- 
ics. Functionalization can be achieved either through the 
adsorption of foreign atoms at different coverage or sub- 
stitution of foreign atoms for B or N in honeycomb struc- 
ture. We considered a number of foreign atoms, such as 
Sc, Ti, V, Cr, Mn, Fe, Mo, W, Pt, H, C, Si, B, N, O, 
Ca, Cu, Pd, Ni and Zn. Part of these atoms are bound 
with a significant energy and form chemical bonds with 
2D BN. Cu, Fe and Si have binding energy smaller than 
0.25 eV, but Cr, Mn, Mo, W, H, N, Ca and Zn cannot 
bind to 2D BN. Owing to the van der Waals interaction, 
the true binding energies can be 0.1 — 0.2 eV larger than 
those calculated with GGA in the present work. 

High coverage of adatoms corresponding to 6 = 1/8 
leads to dramatic modifications in electronic structure, 
if the related binding energy is significant. Under these 
circumstances, either the wide band gap of 2D BN can be 
reduced or diminished and the system becomes metallic. 
Under certain circumstances, the nonmagnetic 2D BN 
attains magnetic moment. Ni, Pd, and Pt covered 2D BN 
are nonmagnetic semiconductors with band gap relatively 
smaller than that of the parent 2D BN. Sc, Ti forming a 
(2x2) structure on 2D BN are AFM metals, V+2D BN 
is an AFM semiconductor. Remarkably, B+ and C+2D 
BN undergo a (4x4) reconstruction and have a band gap 
smaller than the parent 2D BN. Oxygen covered 2D BN is 
found to be a FM, small band gap semiconductor, which 
may display high spin-polarization under bias voltage. If 
the interaction between adatom and 2D BN is weak, the 
electronic structure can be viewed as the combination 
of electronic band structure of 2D BN monolayer and 
adatom monolayer. Si and Fe are weakly bound to 2D 
BN and change the wide band gap of 2D BN to a half- 
metal. 

At low coverage corresponding to O = 1/32, the large 
distance between adatoms hinders any significant inter- 
action between them. This situation is taken to mimic a 
single, isolated atom adsorption to 2D BN, which gives 
rise to localized states in the wide band gap of 2D BN 
honeycomb structure. This conjecture is confirmed by 
examining the adsorption of single C, O, Sc, Ti and Pt to 
(8x8) supercell resulting in ~ 20 A adatom-adatom dis- 
tance. In this respect, the adsorbed adatoms at O = 1/32 
act as dopants of 2D BN. In a few cases the adsorption 
site and magnetic state undergo a change by going from 
the high, O = 1/8 to low coverage, O = 1/32. 

Not only adatoms, but also substitution of foreign 
atoms for B or N in the honeycomb structure give rise 
to localized impurity states in the gap, which attribute 
interesting electronic properties to the system. In partic- 
ular, the substitution of C for B yields an excess charge 
and gives rise to two donor states near the edge of con- 
duction band, the lower lying spin- up band being full. 
In contrast, the substitution of C for N yields a single- 
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Figure 6. (Color online) Schematic energy level diagram of Al and C substituting B; and C and P substituting N in a (4x4) 
supercell of 2D BN. The light/grey shaded region in the background is the valence and conduction band continuum of 2D BN. 
Zero of energy is set to the Fermi level of bare 2D BN. Spin- up bands are indicated with dark/blue lines, whereas spin-down 
bands are light/green. Thick solid bands (levels) are states localized at the substituted atom and nearest atoms of 2D BN. 
Thin solid bands stand for the delocalized states. The energies of some of the localized states with respect to the valence band 
edge are shown in units of eV. Charge density isosurfaces of some of the relevant states are shown below the corresponding 
band diagram. The isosurface value is taken as 7xl0 -5 electrons/A 3 . 



electron deficiency and gives rise to two acceptor states 
above the top of valence band, the lower lying spin-up 
state being full. In both cases, the unpaired spins of C 
atom give rise to magnetic moment \i = 



In conclusion, while 2D BN is a mechanically stiff and 
nonmagnetic wide band gap semiconductor, its band gap 
can be engineered through adatom decoration. In specific 
cases 2D BN attains magnetic properties and becomes 
metallic. Foreign atoms adsorbing at low coverage or 
exchanging with B or N atoms give rise to donor like or 
acceptor like states in band gap. At the end, the material 
achieves interesting properties. Some of these properties 



can be exploited in future applications. 
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